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Analytic models for hadron-hadron scattering are characterized by simple analytical parametriza- 
tions for the forward amplitudes and the use of dispersion relation techniques to study the total 
cross section atot and the p parameter (the ratio between the real and imaginary parts of the forward 
amplitude). In this paper we investigate simultaneously four aspects related to the application of 
the model to pp and pp scattering, from accelerator to cosmic-ray energies: (1) the effect of dif- 
ferent estimations for atot from cosmic-ray experiments; (2) the differences between individual and 
global (simultaneous) fits to atot and p; (3) the role of the subtraction constant in the dispersion 
relations; (4) the effect of distinct asymptotic inputs from different analytic models. This is done 
by using as a framework the single Pomeron and the maximal Odderon parametrizations for the 
total cross section. Our main conclusions are the following: (1) Despite the small infiuence from 
different cosmic-ray estimations, the results allow us to extract an upper bound for the soft Pomeron 
intercept: 1 -\- e = 1.094; (2) although global fits present good statistical results, in general, this 
procedure constraints the rise of atot', (3) the subtraction constant as a free parameter affects the 
fit results at both low and high energies; (4) independently of the cosmic-ray information used and 
the subtraction constant, global fits with the Odderon parametrization predict that, above y/s ~ 70 
GeV, ppp{s) becomes greater than Ppp(s), and this result is in complete agreement with all the data 
presently available. In particular, we infer ppp = 0.134 ± 0.005 at ^/s = 200 GeV and 0.151 ± 0.007 
at 500 GeV (BNL RHIC energies). A detailed discussion of the procedures used and all the results 
obtained is also presented. 

PACS numbers: 13.85.Dz, 13.85.Lg, 13.85.-t 



I. INTRODUCTION 



High-energy soft processes are presently a topical problem in high-energy physics, mainly because they are essentially 
a nonperturbative phenomenon ^} . Despite important results at the interface between soft and hard physics and recent 
progress in nonperturbative QCD, elastic scattering, the simplest soft process, cannot be described in a pure QCD 
framework, but only in a phenomenological context. 

In this area, a variety of models, such as, for example, Regge, Diffraction, QCD-inspired models, and others 0, 
has survived due to a solid theoretical basis and efficiency in describing physical quantities. What is presently 
well established and accepted concerns some general principles, limits, bounds, and theorems, earlier deduced from 
the Mandelstam representation, potential scattering, and also from axiomatic field theories 0, 3- In this context, 
analyticity, unitarity and crossing play central roles and are also the framework of several models referred to above. 
Among them, the so called analytic models are characterized by parametrizations of the hadronic amplitude through 
general analytic functions that strictly obey the formal principles and theorems. Specifically, the aim is to investigate 
two fundamental physical quantities that characterize the forward elastic scattering at high energies, namely, the total 
cross section atot and the p parameter (the ratio of the real to the imaginary part of the forward elastic scattering 
amplitude). In terms of the scattering amplitude F they may be given by 



lniF(s,t = 0) ReF(s,i = 0) 

O-tot S = , P(S) = TT— -r, (1 

s Im }< [s,t — 0) 

where t is the four-momentum transfer squared and s the center-of-mass energy squared. These expressions, in terms 
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of real and imaginary parts of the amplitude, obviously suggest dispersion relations as a suitable formal framework for 
investigations. For particle-particle and particle-antiparticle interactions, the addition of crossing symmetry extends 
and completes the analytical approach. It is expected that such a general formalism, avoiding details of the interaction 
or dynamics, could be a suitable tool in the search for adequate calculational schemes in nonperturbative QCD. 

The analytic models have a long history and important results have been obtained through both integral relations |^ , 
for example and derivative (analyticity) relations for example 0,0,0]. Recently, much effort 

has been concentrated in the COMPETE Collaboration (Computerized Models, Parameter Evaluation for Theory 
and Experiment), which joined the COMPAS group and others specialists also with outstanding contributions in the 
area. These authors have investigated a large class of analytic models, through several statistical indicators that 
complement the usual and and C.L. criteria. One of the main results is the universality of the B\n^ s/sq increase 
of the total cross section for all the collisions considered 0, . 

Despite all the experience accumulated and the detailed analyses that have been developed, we understand that 
three aspects yet need some investigation and this is the main purpose of this work. These aspects are based on the 
following observations. 

1. Beyond the energy region of the accelerators, experimental information on pp total cross sections exists from 

cosmic-ray experiments. Despite the model dependence involved, the large error bars in the numerical results, 
and also the existence of discrepant values, some of these results are usually displayed as a support for several 
model predictions. However, this set of experimental information is not usually taken into account in explicitly 
quantitative analyses that could provide, for example, bounds for the increase of atot and p (and, consequently, 
the intercept of exchanged trajectories) or for the differences between pp and pp total cross sections (for excep- 
tions, see 0,0). Moreover, as will be discussed in some detail, several results obtained by different authors, 
through different approaches, indicate nearly the same increase of crtotis) and this increase is faster than usually 
believed or accepted. 

2. In general, the fits are performed with the full hadronic amplitude, that is, simultaneous fits to atot and p. 

Although some authors correctly claim that this procedure maximizes the number of data points, it must first 
be recalled that important results have been obtained through fits to only the (Ttotdata; a classical example is 
the approach by Donnachie and Landshoff [IS] . Moreover, as will be discussed, in the bulk of experimental 
data available, trtotand phave different status as physical quantities, since pis estimated either by extrapolations 
through dispersion relations (and fits to CTtotdata), or as a fit parameter to the differential cross section in 
the region of Coulomb-nuclear interference, and this is a delicate problem. Finally, since from the previous 
observation we shall also be interested in the high-energy cosmic-ray results, which concern only utot (and not 
p) , global fits may constrain the possible increase of this quantity. 

3. The connection between atot and p through standard dispersion relations demands one subtraction 13 . Although 

the subtraction constant works as a fit parameter in traditional analysis ^ 2^ 8, 9], it does not appear in the 
approach by the COMPETE Collaboration, since the derivative relation or prescription used does not involve 
subtraction constant. Once the analytical approach demands fits to experimental data and the free parameters 
involved are all correlated, it is expected that the presence or not of the subtraction constant may lead to 
different results. 

The aim of this work is to investigate the above three observations and their inter-connections in a quantitative way. 
To this end, we shall analyze only pp and pp elastic scattering, since for particle and antiparticlc interactions they 
correspond to the highest energy interval with available data and are the only set including the cosmic-ray information 
on total cross sections {pp scattering). We shall use as framework two well known analytic models characterized by 
distinct asymptotic inputs, the Pomeron-Reggeon model by Donnachie and Landshoff IT^] and the maximal Odderon 
model in the form discussed by Kang and Nicolescu jl2| ] . For the connections between real and imaginary parts of the 
forward amplitude we shall use derivative (analyticity) relations with one subtraction. Observation 1 is treated by the 
selection of two different ensembles of data from cosmic-ray results, observation 2 through individual and simultaneous 
fits to (Ttot and p, and observation 3 by treating the subtraction constant as a free fit parameter or assuming that its 
value is zero. 

The manuscript is organized as follows. In Sec. II we present the experimental data to be analyzed and the criteria 
for the selection of two ensembles of data. In Sec. Ill we review the main formulas in the analytical approach, 
including the dispersion relations, the models to be used and some high-energy theorems. The fits and results for all 
the cases considered are presented in Sec. IV and the conclusions and final remarks in Sec. V. 
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II. EXPERIMENTAL INFORMATION AND ENSEMBLES 



In this section we review the experimental information that constitutes our data base (accelerator and cosmic-ray 
regions), together with discussions concerning the differences in the determination of (Tjotand p, as well as the model 
dependences and discrepancies related to all the cosmic-ray information presently available. The criteria for selecting 
two ensembles of data are also presented and justified in detail. 



A. Accelerator data 



In the case of accelerator experiments, data on atot from pp and pp scattering and extracted values for the p 
parameter have been accumulated for a long time. Presently, experimental information extends up to 62.5 GeV and 
1.8 TeV for pp and pp scattering, respectively. The database, analyzed and compiled by the Particle Data Group 
(PDG), has become a standard reference and the corresponding readable files are available 0|. Since recent analysis 
has shown that general fits to these data are stable for ^/s above ~ 9 GeV 0| , we shall use here the sets for energies 
above 10 GeV; in our analysis the statistic and systematic errors are added linearly. 

At this point, we briefly recall some differences in the determination of crtot and p 2], which suggest that they do 
not have the same status as physical quantities; this, in turn, corroborates our motivation to investigate global and 
individual fits separately. 

First, one of the methods that was used to extract (Jtot^^ the ISR did not depend on the pvalue, but only on 
the machine luminosity and the rates of elastic and inelastic interactions. Therefore in these determinations (Xtotand 
pare independent quantities. Other methods demand the determination of pand the elastic scattering rate at small 
momentum transfer, extrapolated to the forward direction (slope parameter). In these cases, the quantities to be 
determined correspond to cr^ot(l -t- p^) (luminosity dependent) or atot{"^ + p^) (luminosity independent) "^l. In both 
cases, since it is known that p < 0.14, atotT^sy be obtained with reasonable accuracy even with a rough estimation 
of the ^parameter. Specifically, some authors use the pvalue extrapolated from fits to atot&^d dispersion relations; in 
the other procedure, pis determined from fits to the differential cross section data in the region of Coulomb-nuclear 
interference. In this case the determination is model dependent and it is interesting to note that the procedure, in the 
last instance, demands knowledge of how the hadronic exchanged object interacts, which is exactly what is looked for. 
This point seems clear in the Donnachie-Landshoff approach, since the authors do not use the pdata as input. We 
understand that all these facts reinforce the differences in the determination of pand atot , putting some limits on the 
interpretation of dtotand pas physical quantities with the same status. Beyond these motivations for discriminating 
between individual and global fits, we add our interest in investigating cosmic-ray information, which concerns only 
(Ttotand not p. 



B. Cosmic ray information 

For pp collisions, the total cross section may also be inferred from cosmic-ray experiments and estimations exist in 
the high-energy interval =6 — 40 TeV. The procedure is model dependent and different analyses lead to different 
results, as briefly reviewed in what follows. 

The extraction of the proton-proton total cross section is based on the determination of the proton-air production 
cross section from analysis of extensive air showers. Detailed reviews on the subtleties involved may be found in Refs. 
[23,1^ 1121. Here we recall only the two main steps, stressing the model dependence involved. 

The first step concerns the determination of the proton-air production cross section, namely, the "inelastic cross 
section in which at least one new hadron is produced in addition to nuclear fragments" |2H . This is obtained by the 
formulas 



prod I , N 2.4 X 10"* / , 2x Aatt 

<-air (mb) = — 27' V"-- (g/cm ) = -7-, 

Ap-air (g/cni ) k 

where \p-air is the interaction length of protons in the atmosphere, \att is the shower attenuation length, and the 
inelasticity coefficient A: is a measure of the dissipation of energy through the shower. Aatt is an experimental quantity 
determined through the Xmax attenuation method or the zenith angle attenuation technique |23|. On the other hand, 
the coefficient fc is model dependent and obtained through Monte Carlo simulation; roughly k w 1.5 ^ « 1 (fc > 1) 
when going from Fcynman scaling models to strongly scaling violation models in the fragmentation region 21]. 
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In a second step af^^ is obtained from (7^°^^ through the multiple diffraction formalism (MDF) by Glauber and 
Matthiae 0|, and taking into account the different processes and effects in the p-air total cross section: 



^p-air ^p-air ^p-air ^ p-air 

where (JpZ^l^ concerns the quasielastic excitation of the nucleus and Act is the Gribov screening correction due to 
multiple scattering pnl ]. 

Two important inputs at this point are the nucleon distribution function and a relation between the slope parameter 
B and ct^^j, necessary in the parametrization of the scattering amplitude: 



Fpp{s,t) CX CTfPjCXp 



B{s)t 



With the screening correction, the MDF allows the determination of all the above cross sections and from the strong 

prod 

7 

p-air 



correlation between (7^-°fj. — (jfot ~ ^ the pp total cross section may eventually be estimated (for more details, see 



1211123). 

The results presently available from cosmic-ray experiments, in the above energy region, are characterized by 
discrepancies, mainly due to the model dependence of k and B{s). Before reviewing these results, it is important to 
draw attention to two facts. First, as is well known, discrepancies also characterize some CTtotdata from accelerator 
experiments at the highest energies, for example, at 541 — 546 GeV and mainly at 1.8 TeV. Second, when performing 
quantitative analyses, it is fundamental to select as complete a set of information as possible, or at least those obtained 
in similar circumstances or bases. In the case of cosmic-ray information, it is important to stress that there is no 
experimental determination of af^j, since all the results are model dependent. Therefore they cannot be distinguished 
in terms of an "experimental" status and all the information available must be considered. 

Despite the model dependence involved, we can classify the complete set of cosmic-ray information available accord- 
ing to the inputs or procedures used and, simultaneously, by the value of the total cross section extracted, as explained 
in what follows. 

On one hand the results usually quoted in the literature concern those obtained by the Fly's Eye Collaboration 
in 1984 and Akeno (AGASA) Collaboration in 1993 27]. In the first case, the authors used k = 1.6 (Feynman 
scaling model), a relation of proportionality between B and af^t, namely, the geometrical scaling model [23, and a 
Gaussian profile distribution for the nucleus. At ^/s = 30 TeV, from the value CTp™j^ = 540 ± 50 mb, they extracted 
afot — 120 ± 15 mb 26]. The Akcno Collaboration used k = 1.5 (Feynman scaling model) and the Durand and Pi 
approach to extract the proton-proton cross sections in the region ^/s ^6 — 20 TeV. We shall return to these 
results later on. 

On the other hand, in 1987 Gaisser, Sukhatme, and Yodh (GSY) making use of the Fly's Eye result for CTp™,^, as an 
estimate of the lowest allowed values for CTf^^, and the Chou-Yang prescription for B{s) 31] obtained ct*°* = 175 j^j? 
mb at -y/s = 40 TeV 22]. In 1993, Nikolaev claimed that the Akeno results should be corrected in order to take into 
account the differences between absorption and inelastic cross sections, leading to an increase of the results by « 30 
mb 1^3 • The analysis by Nikolaev was also motivated by previous results from a QCD model of the Pomeron [33l |. 

At this point some fundamental facts concerning all these results must be stressed. First, as reviewed above, the 
Fly's Eye and Akeno results for af^^ are as model dependent as those by Nikolaev and GSY and therefore cannot be 
considered as experimental results; what can and must be discussed is the class of model used in each case. In this 
sense, we first note that since the Fly's Eye Collaboration used the geometrical scaling hypothesis, which is violated 
even at the collider energy, their result is probably wrong. Second, in 1990, Durand and Pi asserted in Ref. [s^ 
that their results published in 1988 '29|, and used by the Akeno Collaboration [23, should be disregarded due to a 
wrong approximation concerning fluctuations. The new results presented in [30j (Fig. 11 in that paper) introduced 
significant changes in the scenario. For example, for ap-air ~ 550 mb, the Akeno result with the method by Durand 
and Pi 29] was ct^^j = 124 mb and from Fig. 11 in Ref. 's^ the value extracted is af^f. = 137 mb, that is, an increase 
of w 11%; to the Fly's Eye published point 26], the corresponding value is a^^^ = 135 mb i3Qj, that is, an increase of 
~ 13%. However, to our knowledge no other publication by Durand and Pi on the subject appeared in the literature. 
Finally, it should be stressed that the Chou-Yang prescription, used by Gaisser, Sukhatme, and Yodh, is somewhat 
model independent and, also to our knowledge, the results by these authors together with those by Nikolaev have 
never been criticized, despite the intrinsic model dependence involved. In conclusion, we understand that all these 
facts, not usually discussed in the literature, suggest an increase of the total cross section faster than indicated by the 
Akeno and Fly's Eye results and that this indication has a reasonable basis. 

It should also be recalled that, more recently. Block, Halzen, and Stanev (BHS) obtained estimations for a^^^ 
from ap-air [l^l, through a QCD-inspired model I35I] and the MDF. In the analysis, the inelasticity coefficient k is 
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TABLE I: Estimations for o-fJ((s) (in mb) from cosmic-ray data [l9l |. 



Elab 




Akeno (AGASA) 


Fly's Eye 


Nikolaev 


GSY 


BHS 


eV) 


(TeV) 


[27] 


[26J 


[32j 


[22] 


[37] 


2.02 


6.2 


93±14 


- 


120±15 


- 


91± 15 


3.52 


8.1 


101±16 


- 


130±18 


- 


100± 18 


(\ 1 1 


in fi 


1 1 7+ 1 8 




1 "14+ 1 7 




1 1 S+ 1 7 


10.63 


14.0 


104±26 




135±29 




103± 29 


18.47 


18.4 


100±27 




129±30 




99± 30 


32.09 


24.3 


124±34 




162±38 




124± 37 


47.96 


30.0 




120±15 






103± 22 


85.26 


40.0 








175±34 





considered as a free parameter, determined through a global fit, including pp and pp accelerator data and the published 
results on the p-air cross section |34|. Specifically, k is determined in a way that their predictions for crf^j = crf^j 
(asymptotically) matches ap-air through the MDF. The resulting value k k, 1.35 seems to be in accord with those 
obtained by combinations of different simulation programs, namely, k ~ 1.15 — 1.30 36]. The extracted af^^{s) at the 
cosmic-ray energies shows agreement with the Akeno results and is about 17 mb below the Fly's Eye value at 30 TeV. 

The numbers associated with all the above cosmic-ray estimations of crf^j are displayed in Table 1. Despite the 
importance and originality of the approach by BHS, we shall not include it in our analysis, because they are in 
agreement with the Akeno results (Table 1) that will be used. However we will treat this case in further work. With 
the exception of the BHS results, all the other cosmic-ray estimates discussed above together with the accelerator 
results are displayed in Fig. 1. 
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FIG. 1: Total cross sections {pp and pp above -^s = 10 GeV): accelerator data and cosmic-ray information available (Table I). 
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C. Ensembles 

From Fig. 1 we see that, despite the large error bars in the cosmic-ray region, we can identify two distinct set 
of estimations: one represented by the results of the Fly's Eye Collaboration together with those by the Akeno 
Collaboration; the other by the results of Gaisser, Sukhatme, and Yodh with those by Nikolaev, which follow the 
higher estimates by Durand and Pi [s^ ■ Taken separately these two sets suggest different scenarios for the increase 
of the total cross section, as already claimed before 0, Is^ . 

Based on these considerations, we shall investigate the behavior of the total cross section by taking into account 
the discrepancies that characterize the cosmic ray information and, to this end, we shall consider two ensembles of 
data and experimental information with the following notation: 

• Ensemble I: pp accelerator data and pp accelerator data + Akeno + Fly's Eye; 

• Ensemble II: pp accelerator data and pp accelerator data + Nikolaev + GSY. 

To some extent, ensemble I represents a kind of high-energy standard picture and ensemble II a nonstandard one. 
However, the consistence among the results by Nikolaev and GSY must be stressed , in addition to their agreement 
with the last results by Durand and Pi (Fig. 11 in Ref. i3^). We add also that, from the discussion in Sec. II B, 
both ensembles seem equally probable. 

III. ANALYTICAL APPROACH 

In this section we first review the essential formulas in the derivative dispersion approach and recall some funda- 
mental theorems and high-energy bounds. The parametrizations that characterize the analytical models to be used 
in the next section are also presented. 

A. Analyticity relations 

For pp and pp scattering, analyticity and crossing symmetry allow us to connect <Jtot{s) and p{s) through two 
compact and symmetric formulas: 



p^PaZ{^^) = E{a+) + 0{a^), 



(2) 



E{a+)-0{a-), 



(3) 



where 



CT±(s) 



-I- a'PP 
"tot 



(4) 



and £'((t+),0(ct_) are analytic transforms relating the real and imaginary parts of crossing even and odd amplitudes, 
respectively. These analyticity relations are usually expressed in an integral form and in the case of the forward 
direction the standard once subtracted integral dispersion relations may be expressed by 0, Q 



E,^t{<J+) = — + — ds' 





1 


/ ds' 


g/2 _ g2 


' So 



a+{s') = 



Re F+{s) 



(5) 



2 

TT 





■ s' ' 






gl2 _ g2 





Re F_(s) 



(6) 



where K is the subtraction constant. Since we shall be interested in the high-energy region, we used in the above 
formulas the cm. energy, instead of the laboratory energy and momentum 

At sufficiently high energies integral analyticity relations may be replaced by derivative forms, usually called ana- 
lyticity relations, which are more useful for practical calculations. In these formulas, differentiation with respect to 
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the logarithm of the energy occurs in the argument of a trigonometric operator expressed b y it s series 0, ■ 
Here we shall use the standard form deduced by Bronzan, Kane, and Sukhatme (see also lldls^ ). obtained from 
the integral form in the high-energy limit: 



, , K 

Eder ) = 'r tan 

s 



2 dins 



a+{s) 



Re F+(s) 



Oder{cr-) = tan 













dins/ 



a_(s) 



Re F^{s) 



Operationally these transforms may be evaluated through the expansions 



K 

Eder{<^+j 

S 



Oderi<y-) 



n d 1 
2 dins 3 

d 



2 dins 



dins 

1 - 



TT d 

' 2 dh^ 
1 /tt d 



3 V 2 dln^ 



2 /vr d 
^ 5 V2 dins 

(T_(s)|' dins 
1 /tt d 



45 V 2 dln^ 



(7) 



cr_ (s) > din s. 



(8) 



This completes the analytical approach: with an input paramctrization for (Ttot{s), Eqs. (l)-(4) allow, in principle, 
the determination of p{s), by means of either the integral forms Eqs. (5)-(6) or the derivative (analyticity) ones Eqs. 
(7)-(8). 



B. Theorems at high energies 

For future reference, we briefly recall here two rigorous high-e nerg y results for the asymptotic behavior of the cross 



sections. The Froissart-Martin bound states that as s ^ oo 



;n-energji 

Mm 



cTtot < C'ln^ s, 

and, according to the generalized form of the Pomeranchuk theorem [3ll42|. if the Froissart-Martin bound is reached 
the difference between pp and pp total cross sections goes as 

pp , pp 

\„ — rrPP rrPP <r „ "*ot ^ "tot 

m s 

which means that the difference may increase as Ins. From Eq. (4), this difference is given by the crossing odd 
component 



Aatot(s) - -2a_(s), 

and therefore, if the Froissart-Martin bound is saturated, a rigorous result is that, asymptotically, the difference 
Acrfot(s) is controlled by the odd component and the maximum contribution is given by ct™°^(s) = Ins. This 
corresponds to one of the variants of the Odderon picture |l2l |4^ and the increase as In s is the maximum Odderon 
hypothesis. If the odd contribution in the imaginary part of the amplitude is not present at the highest energies, then 
Adfot — 0. The possible effects in the real part will be discussed in what follows. 



C. Analytic models 

In the formulas that follows we denote s/s' with s' = 1 GeV^ by s. We shall consider two different parametrizations 
for the total cross section, the main difference being the asymptotic limits, which allow the dominance of an even or 
odd amplitude. 
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1. Donnachie-Landshoff 
The Donnachie-Landshoff (DL) parametrizations for the total cross sections are expressed by 



(9) 



where, originally, X = 21.7 mb, Y = 56.08 mb, Z = 98.39 mb, e = 0.0808, and 77 = 0.4525. We observe that this 
model predicts that the difference between the two cross sections is given by 

Act = aZ{s) - aZ{s) = {Z - Y)s-i (asymptotically). 

Through the formalism described in Sec. Ill A, substitution of the parametrizations (9) into Eq. (4) and then in 
Eqs. (5)-(6) with So = or Eqs. (7)-(8) gives 



K 






h 


Xtan(^); 


s' - 


s 







0(cr- 



and from Eqs. (2)-(3) we obtain the analytical expressions for p{s) 
1 











Xtan(^); 


















Xtan(^); 











(^cotf^)-iI±^tan(^ 
2 V 2 / 2 V 2 



(^COt(^)-(I±^tanr^^ 



■ tan , „ , 

2 V 2 y 



Since in this model Act — > asymptotically, for aZ — <^tot = '^tot we have 



(10) 



(11) 



Ap = pPP - pPP ^ -^AZ Y) cot (^) 



as s 



00. 



8. Kang-Nicolescu 

The parametrizations for the total cross sections used by Kang and Nicolescu (KN) under the hypothesis of the 
Odderon j4^] are expressed by 

aZ{s) = Ai+Bi\ns^k\n^ s, ct^^ (s) = + ^2 Ins + fcln^ s + 2i?s-i/^ (12) 

Differently from the previous case, this model predicts that the difference between the two cross sections is given 

by 



Act = aZ{s) - aZ{s) = {A^ - A,) + {B^ - Bi) Ins + 2Rs-''^ 
— > Ay4 + Ai?lns (asymptotically), 



(13) 



so that, if A A 7^ and/or AS 7^ 0, the total cross section difference may increase and aZt ma-Y even become greater 
than CT^ot, depending on the values and signs of AA and Ai?, which is formally in agreement with the theorems of 
Sec. Ill B. 

With a similar procedure as in the previous model the use of the analyticity relations (7)-(8) leads to 
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and from Eqs. (2)-(3) 



^ i {f - ? (^) - - ^) '»•' - (^) ■} ■ 

In this case, if A A and Ai? are sufficiently small, so that we may replace crf^j « cr^gt = '^tot('S), then, asymptotically, 

Ap = pPP - pPP \AA In s + AB In^ s] . (16) 

This means that, depending on the fit results, there may be a change of sign in Ap, with p^P becoming greater than 
pPP at some finite energy. 

All these possibilities for a change of sign in the differences between pp and pp total cross sections and/or the p 
parameters are based on the concept of the Odderon the case of a change of sign in Aat^s) was early discussed 
by Bernard, Gauron, and Nicolescu 44] and that in Ap by Gauron, Nicolescu, and Leader M^. They are associated 
with the condition that the maximal odderon dominates the imaginary or real part of the amplitude, respectively. 



IV. FITTING OF THE DATA AND RESULTS 



In order to investigate all the points raised in Sec. I and their possible interconnections, we perform 16 different 
fits through the program CERN-MINUIT. In these fits we use both ensembles I and II defined in Sec. II and both 
the DL and KN models described in Sec. III. For each of these four possibilities we perform global and individual fits 
to cTtotand pand, in each case, we either consider the subtraction constant K as a free fit parameter or assume K = 0. 

In the individual approach we first fit only the total cross section data and then extract the corresponding p{s) in 
the case oi K = 0, or, with the results for atotis), we fit only the p data with i^T as a free fit parameter. 

The numerical results of the fits and statistical information are all displayed in Tables II and III for the DL and 
KN models, respectively. The corresponding curves together with the experimental data are shown in Figs. 2, 3, and 
4 for the DL model and Figs. 5, 6, and 7 for the KN model. 
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TABLE II: Results with the Donnachie-Landshoff model. 



r It: 




Individual 






Global 




Quantity: 


fftot 


P 




ftot 


and p 




Ensemble: 


I 


II 


T TT 
i il 


T 
1 


II 


I 


TT 
11 


No. of points 


102 


102 


63 63 


165 


165 


165 


165 




0.76 


0.96 


ICC 1 O A 

i.oo i.o4 


i.uy 


1.24 


0.84 


n no 
U.yo 




20.0 ±0.7 


19.3 ±0.7 




oi 1 n /i 
zi.D ± U.4 


01 o _!_ n /I 
Zl.Z zb U.4 


01 /I _!_ n /I 
zi.4 dz U.4 


oi 1 _!_ n /I 
zi.l dz U.4 


Y (mb) 


48 ±5 


46 ±5 




51 ±3 


51 ±3 


67 ±6 


67 ± 5 


Z (mb) 


74± 10 


70 ±9 




85 ±5 


84 ±5 


114 ± 11 


112 ± 11 


V 


0.37 ±0.03 


0.35 ± 0.03 




0.43 ±0.02 


0.42 ± 0.02 


0.48 ± 0.02 


0.47 ± 0.02 


e 


0.088 ± 0.003 


0.091 ± 0.003 




0.081 ± 0.002 


0.083 ± 0.002 


0.083 ± 0.002 


0.084 ± 0.002 


K 






235 ± 32 245 ± 32 








306 ± 54 


307 ± 52 


Figure: 


2(a) and 2(b) 


2(a) and 2(c) 




3 




4 




FIG. 2: Fits to pp and pp total cross section data from ensembles I (dotted curves for pp and dashed for pp) and II (solid curves 
for pp and dot-dashed for pp), through the DL parametrization (a) and the corresponding predictions for p(s) with K = (b) 
and K as free fit parameter (c) — columns 2-5 in Table II. 
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FIG. 4: Simultaneous fits to atot{s) and p{s) through the DL parametrization with K as free parameter and ensembles I (dotted 
curves for pp and dashed for pp) and II (solid curves for pp and dot-dashed for pp) — columns 8 and 9 in Table II. 
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TABLE III: Results with the Kang-Nicolescu model. 



Fit: 




Individual 






Global 




Quantity: 




<ytot 




P 




ftot 


and p 




Ensemble: 


I 


II 


I 


II 


I 


II 


I 


II 


No. of points 


102 


102 


63 


63 


165 


165 


165 


165 


xVdof 


0.78 


0.91 


1.05 


35.2 


0.78 


0.87 


0.77 


0.86 


Ai (mb) 


44±1 


47±1 




_ 


44.7±0.6 


45.2±0.6 


44.4±0.6 


45.0±0.7 


A2 (mb) 


45±3 


52±3 






44.7±0.7 


45.2±0.7 


44.5±0.7 


45.1±0.7 


Bi (mb) 


-2.9±0.3 


-3.6±0.4 






-3.0±0.2 


-3.1±0.2 


-2.9±0.2 


-3.1±0.2 


B2{mb) 


-2.9±0.6 


-4.2±0.6 






-2.9±0.2 


-3.1±0.2 


-2.9±0.2 


-3.1±0.2 


k (mb) 


0.33±0.03 


0.39±0.03 






0.34±0.01 


0.35±0.01 


0.33±0.01 


0.35±0.01 


R (mb) 


24±7 


12±7 






25.9±0.8 


25.9±0.8 


25.3±0.9 


25.4±0.9 


K 






-198±32 


-1369±32 








-72±46 


-63±46 


Figure: 


5(a) and 5(b) 


5(a) and 5(c) 




6 




7 




FIG. 5: Fits to pp and pp total cross section data from ensembles I (dotted curves for pp and dashed for pp) and II (solid curves 
for pp and dot-dashed for pp) , through the KN parametrization (a) and the corresponding predictions for p{s) with K = (b) 
and K as free fit parameter (c) — columns 2 and 3 in Table III. 



13 




FIG. 7: Simultaneous fits to crtot{s) and p(s) through the KN parametrization with K as free parameter and ensembles I (dotted 
curves for pp and dashed for pp) and II (solid curves for pp and dot-dashed for pp) — columns 8 and 9 in Table III. 
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V. CONCLUSIONS AND FINAL REMARKS 

We shall focus our conclusions and present some discussion on the following five points. 

Ensembles I and II. The figures corresponding to fits to the total cross section data show that, in general, the 
results with ensembles I and II do not differ substantially, except in the case of individual fits with the KN model 
[Fig. 5. (a)]. Also, the cosmic-ray information in ensemble II (Nikolaev and GSY) is not described in all the cases. 
This is a consequence of the small number of points and the large error bars in comparison with the accelerator 
data, as well as the choices for the parametrizations (the models used). However, we must stress that, despite being 
a nonstandard result, the cosmic-ray information in ensemble II has a reasonable basis, as discussed in Sec. II B. 
Certainly, more precise data are necessary for a truly conclusive result, but, at present, it seems reasonable and 
interesting to investigate the semiquantitative consequences of this nonstandard possibility. 

Based on these ideas, in what follows, we shall consider the results with both ensembles. However, it must also be 
noted that some of the conclusions that follow are independent of the ensemble used, such as those concerning global 
vs individual fits and the role of the subtraction constant. 

Bounds for the pomeron intercept. With the DL model and the cosmic-ray information in ensembles I and II, we 
may infer some novel limit values for the soft Pomeron intercept a]p(0) = 1 -I- e. From Table II the highest (lowest) e 
value was obtained with ensemble II (I), in the case of individual fits to crtot (global fits to crtotand p), namely, 

Euppor = 0.094 and eiowcr = 0.079. 

Odderon. In the case of individual fits with the KN model. Table HI and Fig. 5 show that with ensemble II the 
model predicts a crossing in <Ttot{s), so that crf^j becomes greater than crf^j above ^/s « 50 GeV: = 5 ± 3 and 
AB ~ (Sec. II C). However, as shown in Fig. 5, in this case the p(s) data are not described {K = or X as a free 
parameter) . 

On the other hand. Table III shows that, in the case of global fits with both ensembles I and II, statistically, = 
and AB = 0, so that Aatot — <^tot ~ '^tot — 0. However, from Figs. 6 and 7, we see that, in both case, K = and 
if as a free parameter, this model predicts a crossing in the p(s) behavior, with ppp becoming greater than ppp at 
i/s R:! 70 — 80 GeV (for K = 0) and ^ ^ 60 — 70 GeV (for K a free parameter), a result in agreement with the 
early fits by Gauron, Nicolescu, and Leader Most important, these results are also in complete agreement with 
all the experimental data presently available on crtotand pand are independent of the ensemble used. This is certainly 
an interesting and important prediction that will be verified at the BNL Relativistic Heavy Ion Collider (RHIC) and 
the CERN Large Hadron Colhder (LHC) 

Individual and global fits. From Tables II and HI, with the exception of the individual fit to ensemble II with the 
KN model and K as a free parameter, the statistical information does not indicate a preference between individual 
or global fits. 

On the other hand, global fits clearly constrain the possible increase of the total cross section. For example, in the 
case of the DL model, with both ensembles I and II and K — 0, we obtained 

0.081 (reduction « 9%) 
0.083 (reduction w 10%) 

As discussed previously, fifotand pdo not have the same status as physical quantities, since the extracted pvalue 
is model dependent. Therefore, in principle, we understand that global fits underestimate the possible rise of crtothi 
the asymptotic region. In our analysis the effects of the global and individual fits depend also on the subtraction 
constant, as discussed in what follows. 

Subtraction constant. As commented before, if K is taken as a free parameter, the fit procedure demands that it is 
correlated with all the other parameters of the model involved. Therefore, in principle, it is expected to have effects 
in both the low- and high-energy regions. Let us discuss our results through individual and global fits. 

In the case of individual fits with the DL model. Fig. 2 shows that the results for p{s) with K — and K free 
are nearly the same above y/s « 40 GeV and that below this energy the predictions are quite different. We shall not 
discuss the corresponding results with the KN model, since the data are not described (Fig. 5). 

In the case of global fits with the DL model, beyond the differences in p(s), here below ^/s w 40 GeV, the asymptotic 
values of atot{s) are different, namely, from Table II, 



el ,. « 0.088 ^ e^. , 
mdi simul 



« 0.091 ^ , 
mdi simul 
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= 0.081 ± 0.002 and e 



I 

K free 



0.083 ± 0.02 (increase w 2.5%), 



K=0 



= 0.083 ± 0.002 and e 



n 

K free 



0.084 ± 0.02 (increase w 1.2%). 



In particular, with ensemble I, for pp scattering at y/s = 14 TeV the fits indicate af^f ^ 101 nib for A' — and 
crfot ~ 104 mb for K as a. free parameter. On the other hand, for global fits, the KN model is not so sensitive 
to the influence of the subtraction constant at least for p{s) above y/s « 20 GeV, as shown in Figs. 6 and 7. 
Therefore, in general, the subtraction constant affects the fit results in both the low- and high-energy regions. Since 
it is mathematically justified in order to control the convergence of the integral (or derivative) dispersion relation, we 
understand that the subtraction constant cannot be neglected in the analytical approach. 
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